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Abstract 

Using index-free notation, we present the diagonal values ) of the first five heat 

kernel coefficients ') associated with a general Laplace-type operator on a com- 

pact Riemannian space without boundary. The fifth coefficient 05(0;, x) appears here 
for the first time. For the special case of a flat space, but with a gauge connection, the 
sixth coefficient is given too. Also provided are the leading terms for any coefficient, 
both in ascending and descending powers of the Yang-Mills and Riemann curvatures, 
to the same order as required for the fourth coefficient. These results are obtained by 
directly solving the relevant recursion relations, working in Fock-Schwinger gauge and 
Riemann normal coordinates. Our procedure is thus noncovariant, but we show that 
for any coefficient the 'gauged' respectively 'curved' version is found from the corre- 
sponding 'non-gauged' respectively 'flat' coefficient by making some simple covariant 
substitutions. These substitutions being understood, the coefficients retain their 'fiat' 
form and size. In this sense the fifth and sixth coefficient have only 26 and 75 terms 
respectively, allowing us to write them down. Using index-free notation also clarifies 
the general structure of the heat kernel coefficients. In particular, in flat space we 
find that from the fifth coefficient onward, certain scalars are absent. This may be 
relevant for the anomalies of quantum field theories in ten or more dimensions. 



1 Introduction 



The heat kernel method has become a ubiquitous tool in both mathematics and 
physics (see for a recent overview). In mathematics it appears e.g. in the study 
of the spectral geometry of a Laplace-type differential operator on a Riemannian 
space and in the proof of index theorems [0, ^. In physics, the euclidean one-loop 
effective action for a given quantum field theory can be expressed in terms of the 
determinant of such a differential operator which in turn can be written in terms 
of the associated heat kernel. The heat kernel therefore appears in many places, 
from quantum gravity |Q to chiral perturbation theory |jl6|. Anomalies can also be 
studied with the heat kernel method (see |16|, 0). Physicists frequently refer to it as 
the Schwinger-DeWitt @, |] or proper-time method. Exact expressions for the heat 
kernel exist only for special spaces. In the general case one may use its asymptotic 
expansion in the proper-timeQ. The coefficients in this expansion are the so-called 
heat kernel coefficients (see sect 2 for a precise definition). 

Several methods have been developed to find the heat kernel coefficients (see 
e.g. [13 1 for a review). DeWitt determined the first two coefficients^ with a 



covariant recursive method. Sakai relied on Riemannian coordinates to find the 
third coefficient in the scalar case (i.e. a single scalar field on a curved space). For 
the general case, this coefficient was found by Gilkey []^ using a noncovariant pseudo- 
differential-operator technique. The integrated and traced fourth and fifth coefficients 
for an arbitrary field theory in fiat space were found in Q through the evaluation of 
a noncovariant Feynman graph. Avramidi presented a new covariant nonrecursive 
procedure and found the fourth coefficient for the general case (for the scalar case see 
also |jlO[)- More recently, string- inspired world line path integral methods have been 



used [0 to determine the functional trace of the first eight heat kernel coefficients for 
the case of a matrix potential in fiat space without gauge connection (see also |T2|). 



In this paper the explicit diagonal value of the fifth heat kernel coefficient in the 
general case is presented for the first time. In physics this coefficient is of impor- 
tance e.g. in analysing the short distance behavior and anomalies of ten-dimensional 
quantum field theories (see [0). However, the number of terms in the higher heat 
kernel coefficients grows rapidly, leading one to expect more than a thousand terms 
for the fifth coefficient^. This would seem to preclude writing down this coefficient 
in an intelligable form. Indeed, a computer would appear to be an essential piece 
of equipment in determining and storing the fifth coefficient. Contrary to these ex- 
pectations, we will show here that with a suitable index-free notation one obtains a 
compact expression for this coefficient, containing only 26 terms. 

Using standard matrix notation for the field indices, the heat kernel coefficients are 



^An alternative is the so-called covariant perturbation theory It provides a (par- 
tial) summation of the Schwinger-DeWitt series and can account for nonlocal effects. 

^We frequently abbreviate 'diagonal value of a heat kernel coefficient' to 'coefficient'. 

•^Even restricting to the a single scalar field on a curved space, the j-th coefficient 
already contains 1, 4, 17, 92, 668 i^-'-terms for j ~ 1, 2, 3, 4, 5 (see appendix A of [fsl). 
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scalars and this suggests that it may be possible to write them in a form which is free 
of spacetime indices as well. Since the heat kernel coefficients are covariant, we may 
determine them in a special gauge and adapted coordinates. Following earlier authors, 
we select the Fock-Schwinger gauge and Riemann normal coordinates. Examination of 
the relevant recursion relations then shows that only certain maximally symmetrized 
multiple (covariant) derivatives of the matrix potential and Yang-Mills or Riemann 
curvature tensors appear. The heat kernel coefficients being formal scalars, it turns 
out that we need to keep track of the rank of these tensors only (note the similarity 
with totally antisymmetric tensors, i.e. differential forms). This provides the basis 
for our index-free notation. Using instead a fully covariant method, Avramidi 
has arrived at similar conclusions by expanding in so-called covariant Taylor series, 
obtaining in this way the first four coefficients. Our non-covariant procedure seems to 
be no less efficient and, by maintaining index-free notation and manifest hermiticity, 
yields even compacter answers for the heat kernel coefficients. Thus, in the special 
case of a flat space with a gauge connection, we can also present the answer for the 
sixth coefficient. It contains only 75 terms. 

Using index-free notation also permits us to investigate the general structure of 
the heat kernel coefficients. Without actually solving the recursion relations, we can 
show that certain Lorentz scalars are absent from all coefficients. In flat space, this 
happens for the flrst time in the flfth coefficient. This may therefore be relevant for 
the anomalies of quantum field theories in ten or more dimensions. 

As we already mentioned above, a brute force approach with a computer algebra 
program would have produced an unwieldy result. Our results were obtained without 
the aid of a computer. However, it is relatively easy to program our index-free method 
and we used FORM IBTI and Mathematica to run some checks. 



An outline of this paper is as follows. In sect 2 we recall the main features of the 
heat kernel method alias the Schwinger-DeWitt formalism. In sect 3 we introduce 
our index-free notation and use it to determine flrst the heat kernel coefficients in flat 
space without gauge connection, expanding in powers of either derivatives (sect 3.1) 
or of the matrix potential (sect 3.2). In sect 3.1 we also use our index-free notation 
to prove that certain scalars are absent from all heat kernel coefficients. In sect 4 
we show that the corresponding heat kernel coefficients with a gauge connection can 
be obtained from a simple covariantization process. This involves not only replacing 
partial derivatives by covariant ones but also adding new field strength dependent 
terms. The latter kind of terms are shown to arise only as shifts in the potential and 
in its covariant derivatives. These shifts being understood, the heat kernel coefficients 
do not change their form upon 'turning on the gauge field'. In sect 5 we employ 
Riemann normal coordinates to generalize to a curved space. In particular, we present 
in subsect 5.1 an explicit expansion of the vielbein to all orders in these coordinates 
(such an expansion, usually given for the metric, but only to some flnite order in the 
normal coordinates can be found in many places, e.g. ^). We give a similar result 
in subsect 5.2 for the gauge connection to all orders in normal coordinates. Based 
on this, we flnd in subsect 5.3 the explicit form for any heat kernel coefficient up to 
and including terms of fourth order in the Yang- Mills and Riemann curvatures. The 
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curved coefficients can be obtained from the corresponding gauged but flat coefficients 
via further simple covariant substitutions. To complete the fifth coefficient, we need 
to find the few terms of fifth order in the curvatures. This we do in subsect 5.4 by 
specializing to a locally symmetric space. In sect 6 we present the explicit answers for 
the ffist five coefficients. We indicate how to return to more conventional notation and 
compare with earlier results. Our conclusions are given in sect 7. Several appendices 
follow (in particular, the sixth coefficient in flat space is given in appendix C). 



2 Schwinger-DeWitt formalism 

Consider a set of fields (piix), i = l...n, defined over a compact rf-dimensional 
Riemannian manifoldQ with coordinates x'^, ^ = 1 . . .d and metric g^^ (see appendix 
A for our notation and conventions). The fields are acted upon by a Laplace-type 
wave operator A 

A = -V^-X , = g^'^Vf.V^ (2.1) 

Here the covariant derivative V includes connection terms as needed for the fields (pi. 
X is a hermitian n x n matrix potential (we suppress the field or 'bundle' indices). 
The wave operator A is hermitian with respect to the inner product 



(0,^) = J d'x^ (l>*ilj (2.2) 

For most bosonic gauge field theories of interest one can achieve a wave operator of 
Laplace-type as in ( |2.1| ) by a suitable gauge choice. For fermionic (gauge) fields one 



squares the wave operator to obtain again ( p.l|) (see |jT9| for wave operators not of 
this form). 

Following Schwinger and DeWitt, we introduce the proper-time parameter r and 
define the heat kernel K associated with A by 

d 

(— + A)K(x,x';t) = , K(x,x';0) = IS(x,x') (2.3) 

OT 

where / is the n x n unit matrix and the bi-scalar 6 function is defined by 

J d'^x^6{x,x')(j){x) = (f){x') (2.4) 

for any scalar field 0. As we mentioned in the introduction, an exact solution for the 
kernel K exists only for special spaces. We instead make DeWitt 's ansatz 



K{x,x';t) = (47rr)-'^/2p(x,x')'/'e-'^("'"')/2-^a,-(a;,x')^ (2.5) 

j=Q J- 



J 



^The fields can be considered to be sections of a smooth vector bundle. The heat 
kernel coeiBcients do not explicitly depend on the dimension or signature of space(time). 
We assume space to have no boundary (for the case with boundary, see e.g. pO|). 
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which is known to be an asymptotic expansion in r [|^. Note our unconventional 
normahzation for the heat kernel coefficients aj , which however agrees with [^j . They 
transform as scalar densities of weight —1/2 at both x and x'. The bi-scalar a is the 
geodetic interval (one half of the distance squared between x and x') and satisfies 

= 2a , [a] = (T{x,x) = (2.6) 

where we use Synge's bracket notation to indicate evaluation on the diagonal. The 
bi-scalar V is the Van Vleck-Morette determinant defined by 

V{x,x') = g-'/' det{-a,,,,) g' (2.7) 

where a prime refers to the point x'. It satisfies 

(2a.^V^ + a.%-rf)I?i/2 = , [V] = l (2.8) 

Inserting (|2.5| ) into (|2.3| ) and using ( p.6|) and ( p.8|) , one finds that the heat kernel 
coefficients must satisfy the following recursion relations for j > 

{cr^V^+j)a, = -jV-'/^AV'/^a,.i , M = / (2.9) 

where it is to be understood that a_i vanishes. Note that whereas a and V depend 
only on the metric, the aj are matrix valued and depend in addition on the detailed 
form of the wave operator. The hermiticity of the wave operator implies 

aj{x,x')^ = aj{x',x) =^ [aj]^ = [aj] (2-10) 

To keep this property manifest we introduce the following notation: for any matrix 
valued function F we define 

{F} = F + F^ if F^ F^ else F (2.11) 

Thus {F} = F when F is selfadjoint. Frequently one is interested in the functional 
trace of the heat kernel coefficients 

bj = Ti ttj = ti J d'^x ^/g aj{x, x) (2-12) 

where tr denotes the matrix trace over field indices only. To determine e.g. chiral 
anomalies one would need the non-traced version. In this paper we will determine 
the non-traced diagonal heat kernel coefficients. 



3 Flat space without gauge connection 

In a flat space without gauge connection the recursion relations (|2.9| ) for the heat 
kernel coefficients become 

{x^d^ + j) a, = J (d^ + X)a,_i , [ao] = I (3.1) 
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Here we have set x' to zero (we will not differentiate at x'\ We are mostly interested in 
the diagonal values of the heat kernel coefficients. However, it is easy to see that this 
in turn requires knowledge of some derivatives of preceding heat kernel coefficients 
on the diagonal. Indeed, taking the diagonal value of ( |3.1|) yields 



K-] = [92a,_i] + X [a,_i] (3.2) 

Thus in particular 

[ai] = X (3.3) 

but for j > 1 we must first find a^.i and d'^aj^i on the diagonal. Applying d"^ to ( p.l|) 
and then going on the diagonal gives 

= ([525^0,-2] + X [92a,_2] + 2X'^ [d.a,,^] + X% [a,_2]) (3.4) 

Setting j = 2 and substituting the result in (p.2|), we find 



M = \d^X + X^ (3.5) 

but for j > 2 we require some new diagonal values of derivatives of aj_2. This 
recursive procedure ends after j steps since the diagonal value of any derivative of ag 
vanishes. Differentiating (|3.1| ) n times and taking the diagonal value yields 



(aj-i,(2),n + (Xaj_i),„) (3.7) 



[d^,, . . . d^^aj] = - [d^^ . . . df^^Aaj^i] (3.6) 

To solve these recursion relations in an effective way, we introduce a short hand 
notation for them. Using comma notation for partial derivatives, but writing only 
the number n of uncontracted derivatives taken, we can abbreviate ( ^.6|) as 

_ 3 
J +n 

Here, the index n stands for all partial derivatives taken on the left hand side of ( p.6|) . 
The sans sen/ symbols serve to emphasize that we are using this short hand notation 
and at the same time imply evaluation on the diagonal. The d"^ on the right hand side 
has been abbreviated to a 2 in parenthesis. Of course one must first distribute the 
n derivatives over the factors of the second term in ( |3.7| ) before taken the diagonal 
value. Doing so yields 

^ 7T^ + £ ( „ ) ^(p ^i-i.'^-p) ) (3-8) 

3 ^ "' p=o \P/ 

where the parenthesis around the p plus n — p indices imply total symmetrization. 
We write the result of replacing n by 2n and contracting all derivatives as 

^i.(2n) = — ( ai-l,(2n+2) + II ^ j ^(p aj_i,2n-p) j (3.9) 
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where the parenthesis now imply not only total symmetrization, but also full contrac- 
tion. Here we introduced the following notation: for any functions F{x), G{x) and 
H{x) we define 

F(2n) = F,{^ll^ll...^in^ln)i^) 
F(fcG2„-fc) = F(^,...;.,(0)G,^,^,..^,„)(0)5^^'^^..5'^--'^- (3.10) 
F(,G,H2„_,_,) = F(,,...^,(0) G,,,^,...,,^M //,^,^,^,...^,„)(0) 6^^^- . . . ^^--'^^ 

with similar expressions in case the 2n derivatives are distributed over yet more 
factors. Unless otherwise noted, the use of parenthesis as on the left hand side 



of (|3.1CI| ) indicates that the 2n derivatives are to be totally symmetrized and fully 



contracted. Note that 

F(2n) = (d'TFiO) (3.11) 

so in this case the symmetrization in ( |3.10| ) is superfluous (it is relevant when a gauge 
connection is present: see ( [4.14|) ). The notation of (|3.11|) is due to Avramidi 



(evaluation at the origin is not implied there). We generalized it here by allowing 
distribution of the contracted and symmetrized derivatives over any number of factors. 
Note that e.g. 

F(iGi) = F^(0)G,^(0) 

F(iG3) = F^(0)G,^,,(0) = F(iGi)(2) (3.12) 
F(2G2) = |f^;.(0)G,,,(0) + |f^,(0)G',^,(0) 

The general rule for such a reduction is given in ( [A.9| ). 



3.1 Expansion in derivatives 

We first discuss some generic features of the derivative expansion, which can be 
understood without solving the recursion relations. Using only dimensional analysis, 
it follows that through second order in derivatives the diagonal values of the heat 
kernel coefficients must have the structure 

[a,] = a,X^ +J2 ajkX'-'X,^^X^-'-' + J2'j2 a,kiX''-'X,,X^-'-'-'X,,X'-' + . . . 

k=l k=l e=i 

(3.13) 

The a's are numerical coefficients which, in order for ( p.l3| ) to be hermitian, must 
satisfy 

aj j-k = ajk , ajek = ajki (3-14) 

These expectations are borne out by our results. In particular, all terms in ( |3.13| ) will 
turn out to have nonvanishing coefficients. At fourth order in derivatives, dimensional 
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analysis allows 11 ways of distributing the derivatives, namely 



,^^iU-^ ,U 5 ,flfl-^ ,UU 5 ,^U-^ (3.15) 

X^^flX^lyX^ly , X^^X^^^X^ly , X^^JyX^^X^J, , X^^X^f^^X^/y 

X ^l_,X ^^X ^yX , X^^X^yX^^X^y , X ^^X ^yX ^yX 



Here and through ( p.l9| ) below, we list only equivalence classes of Lorentz scalars. 



where two scalars are considered equivalent if they become equal upon omission of all 
undifferentiated factors X and/or reversing the order of the factors^. We therefore 
do not write the curly brackets as in ( p.llj ). 



We now claim that two of the Lorentz scalars in (|3.15|) can not appear in the heat 



kernel coefficients. To prove our assertion, we need not solve the recursion relations 
explicitly. It suffices to keep track of the way the derivatives get distributed over the 
various factors as the recursion proceeds. Hence, we omit numerical factors as well as 
any undifferentiated matrix X. We may even drop the ordering labels j respectively 
j — 1 in ( ^.8| ) and thus abbreviate aj_„ to a„, where n is the number of derivatives. 
Writing only Lorentz scalars, we thus have the following chain of substitutions 

a a(2) 

2n 

a(2n) a(2„+2) + X! ^(P^2n-p) , n>l (3.16) 

p=l 

2n—p 

X(pa2n-p) X(pa2n-p){2) + X/ ^{p^q^2n-p-q) , p = 1 . . . 2n 

q=l 

etc. We truncate this hierarchy at level 2N in derivatives, i.e. we drop any term 
with more than 2N derivatives (note that for the j-th coefficient, N < j — 1). To be 
definite, consider the case N = 2 where the first few steps are 

a 3(2) 

a(4) + X(iai) + X(2)a (3.17) 
X(4) + X(3ai) + X(2a2) + X(ia3) + X(iXi)a + X(2)a 



We usedQthat X(iai)(2) = X(ia3), see (|3.12|) Iterating once more, we see that at fourth 



order in derivatives a generic heat kernel coefficient can contain only the following 9 
equivalence classes of Lorentz scalars 

X(4) 



^By working with real fields only, we can always arrange X to be a real symmetric 



matrix. The same is true for its derivatives and {-F} as in (2.11) then means that, unless 
_F is a palindrome, we must add to F its transpose. 

Note that in general we have X(pa2n-p)(2) 7^ X(pa2„_p+2)- Thus, before being able to 
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iterate again, we must 'remove the box'. See ( 3.30 ) for an example and app. B for the 



general solution. Since we shall not go beyond N — 2 this complication is irrelevant here. 
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X(3Xl) , X(2X2) , X(2) X(2) 

X(iX2Xi) , X(2XiXi) , X(2)X(iXi) (3.18) 
X(iXiXiXi) , X(iXi) X(iXi) 



Comparison with ( p.l5| ) shows that expressions with overlapping sets of symmetrized 



and contracted derivatives, namely 

X(iX(2)Xi) = X.X^^X, (3.19) 
X(iX(iXi)Xi) = X^X.X^X^ or X(iX(iXi)Xi) = X^X.X^X, 

are absent. This holds for all heat kernel coefficients (recall that we dropped any 
undifferentiated matrix X). The absence of the first (second and third) entry in 
(|3.19D can be first observed in 35 (respectively ag). This may hence be relevant to 



the short distance behavior and anomalies of quantum field theories in ten or more 
dimensions. We conclude that at fourth order in derivatives only 9 of the 11 a 'priori 
allowed scalars appear. This is only a small reduction in the number of terms, but 
at higher orders the relative number of such absent Lorentz scalars increases rapidly. 
At sixth order in derivatives we find that dimensional analysis and hermiticity would 
permit 85 different Lorentz scalars similar to those in ( |3.15|) . However, continuing 
( |3.17| ), we discover that only 53 of these scalars can actually appear in the heat kernel 



coefficients. 

We now give our explicit solution for a^- through fourth order in derivatives. To 
have manifest hermiticity, we use here the convention that in an A^-fold sum with 
ki, . . . ,kN as summation variables (below ki, ^2, • • • = k, i, . . . etc, but skip o), the 
label /ctv+i has by definition the value (fcAr)max — kj^ + 1 (so e.g. in the second line 
below, i = j — k). 

a, = X^ 

+ E-^x-^X(2)X-^ 

j-2j-k-l r)ip 

+ -^x'^-iX(iX™-iXi)X^-i 

kjk + + 1) vfc-iv ^,,-1 

+ g 20- + !)(, + 2) 

j-3j-k-2 ip 

+ E E (^X'^-iX(2)X-iX(2)X^-^ 



k=i i=i 3 + 

3A;(A; + 1)£(£ + 1) v^-iv y^-i 

1 r~, ; A AcoA Ao^A 

(j + l)(j + 2) 
2k{k + lW + m+l 



{X'^~^X(3X"^-^Xi)X^-H) 
(j + l)(j + 2) ^ ^) 

j-4j~k-3j-k-e-2 OL^™ 

+ E E E (^{x^-%)X-^X(iX"-iXi)X^-i} 



A;=l i=l m-- 
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6k{k + + n + 1) ym-lY vn-lv 

(j + l)(j+2) ^ (2 11)^ 



Qk{k + m+ + n + 1) _i y^-ly y„-ly 

T ; Z~, ; A A/i A AoA AtiA 



j-5 j-fc-4 j-fc-£-3 j-k-i-m-2 



+ E E E E (■^x*^-iX(iX™-iXi)X^'-iX(iX"-iXi)X^-i 



fc=l f=l m=l n=l 



I 12/i;(A; + m + + n + 1) wm-iy yp-1y vn-iv wi-l 



+ Oid") (3.20) 

With the restriction of at most four derivatives, this result constitutes an exphcit 
solution for all heat kernel coefficients in flat space without a gauge connection. It is 
complete for j < 3 and yields all terms but one in 34. To find the coefficient of the 
'missing' X(6) term in 84 it is best to expand in powers of the matrix potential. This 
will be presented in the next subsection. Also note that the last two terms in (|3.20|) 



do not appear until a^. Taken together, ( |3.20D and ( p. 311 ) yield the complete answer 



for the first six heat kernel coefficients. We refer to sect 6 for the explicit answers for 
the first five coefficients (to obtain the fiat space results, replace each Z in ( |6.1| ) by 
an X and omit all hats and daggers). The result for ag is presented in appendix C. 



3.2 Expansion in the potential 

We return to ( |3.7| ), now paying attention to the order in the matrix potential rather 
than in derivatives. Thus we are looking for an expansion that starts as 

[a,] = a'j{dy-^X + 0{X^) (3.21) 

a'j being some j-dependent numerical factor. We will show that the recursion relations 
yield an expression of the form 

aj = a'jX(^2j-2) + E E "iMX(2fc-£X^)(2j_2fc-4) 

k e 

+ E E E E E '^iMmrapX(2fc-m^L-rf ^n){2i-2fc-2£-6) + ^(X'^) (3.22) 

k e m n p 

where we used our shorthand notation (the exponents 2i — p and p on the second 
and third factor of the last term also count derivatives). This is not manifestly 
hermitian and verifying hermiticity thus provides a strong check on our results for 
the a' coefficients. Furthermore, at third order in X, ( |3.22| ) shows overlapping sets of 
derivatives. The results of the previous subsection imply that it must be possible to 
remove these overlapping terms. Below we show how this can be done. 

For convenience, we repeat our starting point, eq ( ^.7| ) 

) (3.23) 
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The first term on the right hand side, which is of zeroth order in X, can be eliminated 
by substituting for it from the left hand side, i.e. replace j— >j — >n + 2 and 

contract one pair of derivatives 

j — 1 / \ 
ai-i,n,(2) = J ^^^-^ [ ^j-2,n,(4) + (X a^.a) ,„,(2) j (3.24) 

Iterating this yields 

^ j — q ■'^^ (fc'' 

aj> = (Xaj-fc-l),n,(2fc) = Y. fj+k+n\ (Xaj-fc-i),n,(2fc) (3.25) 

fc=0 q=0 J ^ ^ 'J k=0 [ fc+i ) 

Clearly, each term is now at least of first order in X. We prefer to reverse the order 
of summation and write this as 



a 



k=0 



Y ^Ik (Xafc),n,(2i-2fe-2) (3.26) 



where we defined the combinatorical coefficients C by 

f2j-k+n-l\ 
\ j-k J 

Now separate off the k = term and distribute the derivatives over Xa^ to obtain 



Cjk - f2j-k+n-l\ (^•2'') 



3j> = ^Jo Xn,(2j-2) + X! Y „ C'Jfc X(2j'+jj-.pafc,p) (3.28) 

fc=l p=0 \ P / 

with j' = j — — 1. The set of indices /ii . . is labeled n here to indicate that 
the elements of this set are to be included in the indicated symmetrization, but they 
remain uncontracted. By substituting for a^^p from the left hand side we find 



j~l 2j'+n /^f + n\ 

aj> = C'Jo Xn,(2i~2) + XI „ )C'j"fe'^fcoX{2j'+n-pXp)(2A;-2) 

k=l p=0 \ P J 

j-l k-1 2j'+n /^f 

'^YY Y [ ^ C'"fcC'LX(2i'+n-p(Xa£)p)(2fe') (3.29) 



k=2e=i p=o 



P 



with k' = k — This shows explicitly the terms of second order in X. To proceed 

to third order, we take n = and distribute the derivatives over Xa^. In general 
requires two binomial sums, one for each of the sets of derivatives marked p and 2k', 
and leads to overlapping sets of derivatives as in ( p. 221 ). To avoid this and to allow 
us to lump the two sets of derivatives together, we use that for any functions F{x) 
and G{x) 

F{i G'i)(2) = Gs) 

Fi2 G2)(2) = ^F(^2 Gi) + iF(2) ^(4) (3.30) 
11 



etc. Choosing F = X, G = Xai and taking x = , the first identity shows that we 
can trivially lump the derivatives together in computing 85 and the second identity 
shows how to achieve the same for ag. In appendix B, we show how to do this for 
arbitrary values of j. Restricting for simplicity here to j < 5, we obtain as our final 
result 

[j/2]j-2k 

= TiZHY X(2j~2) + H 51 '^ifcC'fco{X(2j'-pXp)(2fc-2)} 
[ j ) k=l p=0 \ P / 

j-lk-1 2j' 2k'+p /2k' +p\ 

+ ^^^ J2 [ „ ][ „ C'jfcC'LCloX(2j'-pX2fc'+p_gXg)(2£_2) 
k=2 e=l p=0 q=0 \ " / \ y / 

+0(X^) (3.31) 

with j' = j — k — 1, k' ~ k — i — 1 and the C-symbols were defined in (|3.27| ). Note 
that we rewrote the second order terms in manifestly hermitian form, the range of the 
double sum having been correspondingly restricted (we use [n] to denote the integer 
part of n) . Expression ( |3.31[ ) is one of the main results of this paper. 



4 Flat space with gauge connection 

We will show that the heat kernel coefficients in the presence of a gauge connection 
can be obtained from their 'trivial' counterparts without such a connection by making 
simple covariant substitutions of the kind 

df,^... d^^X V(^i . . . V^.)X + F-dependent terms (4.1) 

Here partial derivatives are turned into totally symmetrized covariant derivatives and 
in general there are additional field strength dependent terms. 

Denote the nonabelian vector connection by A. The associated covariant deriva- 
tive and field strength are defined by 

= (9^ + , F^. = [V^, v.] (4.2) 

We take V and thus F to be antihermitian. The heat kernel coefficients satisfy 

{x^V^ + ]) aj = J (V^ + X) a,_i , [ao] = I (4.3) 

To solve these recursion relations, we find it convenient to work in Fock-Schwinger 
gauge 

x^A^{x) = (4.4) 

This is equivalent to the requirement that all partial derivatives of the gauge connec- 
tion vanish upon total symmetrization at the origin, i.e. for j > 1 

Ami,...m.)(0) = (4.5) 
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In particular, the gauge field vanishes at the origin and for any n > 

(9=^)"9-A(0) = (4.6) 

In Fock-Schwinger gauge the recursion relations simplify to (c/eq ( |3.1|) ) 

{x^'d, + j) a, = J {d^ + X) a,_i , [ao] = / (4.7) 

where we defined the differential operator 

X = X + A^^^ + 2A^'d^ + A^'A^ (4.8) 

Let Z be its non-differential-operator part, i.e. 

Z = X + A\^, + A^A^ (4.9) 

We will soon need a covariant expression for the partial derivatives of Z at the origin. 
It is well known and easily verified that in Fock-Schwinger gauge the partial derivatives 
of the gauge field at this point have covariant values given for j > 1 by 

^V...M.(0) = -^i^(./;....,)(0) (4.10) 

Here we use semicolon notation for Yang-Mills covariant derivatives. Eq ( |4.1CI| ) implies 

= ^i^(/;Mi...M.)(0) , J>1 (4.11) 



If we now define new covariant sans serif symbols Y 



Y/ = -^^(m/;.....)(0) (4.12) 



as 



which are to be treated formally as vectors, then we can abbreviate ( [4.10| , [4.11 
^^,....,,(0) = Y/ , A\^,,...,^{0) = Y(/;,) , J > 1 

^',.....^.(0) = Efilw,-^) , J>2 (4.13) 



k=l 



We further note that at the origin, due to ( |4.5| ), we can immediately covariantize the 
partial derivatives of the matrix potential as follows 

X, = X^,...^,(0) = X(^,...^^.)(0) (4.14) 

Here we used the same symbol Xj as in sect 3 to mean now a totally symmetrized 
j-fold covariant derivative of the matrix potential at the origin. Thus, the desired 
covariant expression for any partial derivative of Z at the origin is 

Z, ^ Z,,,...,^iO) = X, + Y(/^,) + X: Q Y(,Y,„,) (4.15) 
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No implicit contractions occur here except for the scalar product between and 
Yj_fc. Replacing j by 2j and contracting all indices yields 



-(2i) 





{Y(,Y,)} (4.16) 



where i = 2j — k is to be understood. In the second expression we used the nota- 
tion of ( p.ll| ) to obtain a manifestly hermitian result. In sect. 5 we determine the 
generalization of (|4.15| , |4.16| ) to a curved space. In that case we do not find a closed 
expression which holds for all values of j . 



4.1 Expansion in derivatives 

Returning to our short hand notation, we have 

(a,-l,(2),n + (X 



J + n 



(4.17) 



Similar to (|3.16| ), we now have the following chain of substitutions (we again omit 
numerical factors and ordering labels for the heat kernel coefficients, but this time we 
keep undifferentiated matrices X) 



a 

a(2n) 



a(2) + Xa 

a(2n+2) + (Xa)(2n) , ^ > 1 



(4.18) 



etc. In the first line we used that X at the origin equals X. To remove X from the 
second line of ( ^4.18| ) as well, we use that at the origin 



(Xa)(2.) = (Zta) 



(2n) 



n>0 



(4.19) 



Note that the right hand side no longer contains a differential operator. To prove 
this, we use X = Z + 2A^d^ and Z'^ = Z — 2A'^^^ and note that for any function F 
one has (see appendix D for the proof) 



{dyd^{A''F){0) = 



n>0 



(4.20) 



Taking for F a heat kernel coefficient, (|4.19|) follows. Adding one more step to the 
hierarchy ( 4.18 ), we obtain 



3(2) 



+ Xa 



2n 



3(2n) 



■ a(2n+2) 

+ E 

'Z-{p32n-p) 'Z-(p 32n-p)(2) + Z(p^ (X a)2n-p) 



n > 1 



(4.21) 



p = . . .2n 
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Except for p = and p = 2n, ([4.19|) can not be used to eliminate X from the last 
term. Instead we proceed as follows (write r for 2n—p and keep binomial coefficients 
here) 



q=Q yij 



\Ap'^q^r-q) \ \Zi^^2N'q> + iar-q'-l)u 



9=0 q^O \(1 + 
= X! ]'^{p\'^q^r-q) -\ — "^^^ q+l3r-q-l)iy 

= iltViP^^r-q) (4.22) 

In the second line we shifted the summation indexQ q so as to collect terms with the 
same number of derivatives on the heat kernel coefficient. In the third line we can 
use ([4.101) to replace A^+i by Yg+i and thus obtain a covariant answer. The last line 
defines Zq in terms of the previous line. The hatted Z is designed to absorb the new 
field strength dependent term (this term exists even for g = and we will write Z for 
Zq. In case r = too, Zq = X). With this definition of Z we can replace the last line 
of the hierarchy (|4.21| ) by 

2n—p 

'Z-{p32n-p) Z(p^a2„-p)(2) + X! Z(p^^g32n-p-g) , p = . . . 2n (4.23) 

q=0 

Iterating once more does not bring new features^ and we arrive at the following con- 
clusion. 

The diagonal values of the heat kernel coefficients in the presence of a gauge connec- 
tion are obtained from those without this gauge connection by making the following 
covariant substitutions in 2(\ ) 



X(2i) -> Z(2j) 



XqX„X,) ^ Z(;Z„Zfc) (4.24) 
(^jXpXnXk) Z(/ZpZ„Zfc) 



etc, where j, k > 1, n,p, . . . > 0, Zj is defined in ( \4-13i) and the Z act to their right 
as follows 

2 

ZrJ-k) = Z„Zfc) + — kYn+l'^^k-l)u 

n + 1 

^ ^ ^ 2 

ZpZ.„Zfe) = ZpZ„Zfc) H — Yp+/(nZ„_i''Zfc) + A;Z„Zfc_i);,) (4.25) 

n + 1 



^The q' = —1 term is absent because the gauge connection vanishes at the origin, see (4.5) 
^ Other than those remarked upon in footnote 6. 
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Since they are unaffected, we omitted Z(j^ . . . from tlie left of each term in ([4.25|) , the 
dots representing possible further Z. Note that Z„ acts on each one of the k indices 
on Zfc in the same way, hence the factor k. Similarly Zp acts on n plus k indices, etc. 
Thus, upon 'turning on the gauge field', ( |3.20| ) becomes 

i~l up j-2j-k-l r)l.p 

a, = X^ +E-XtX'"'Z(2)X^"1 + ^ ^ ^X'="^Z(,tZ-%)X^-i + ... (4.26) 
k=i3 + ^ k=i e=i J + ^ 

Note that only an X which is 'sandwiched' between X(j and X^.) gets replaced by a 
Z. Formally, the substitutions in ( [4.24 ) do not change the number of terms or their 
numerical coefficients. In this sense the heat kernel coefficients retain their original 
appearance in the gauging process. If we use the curly bracket notation of ( p.ll[ ), we 
should take note that for j k 

{Z(;Z„Zfe)} = Z(;Z„Zfc) + Z(fctZ„Z,-) (4.27) 

with no dagger on Z„ in the second term (this can be shown to be hermitian). Finally, 
repeating the steps of ( p.l7| ), we find that the covariant analogues of ( p.l9| ) are absent. 



4.2 Expansion in the potential 



Upon replacing X in sect. 3.2 by the differential operator X and retracing our steps, 
we find that the heat kernel coefficients for the case with a gauge connection are 
obtained from those in (|3.31|) through the same substitutions as in (|4.24| ). Thus 



\j/2] j-2k 

-(2i-2) + H H 
k=l p=0 



"If 



C^jk^kO {Z(2j'-p'Zp)(2fc-2)} 



j-1 k-1 2j' 2k'+p 

+ EEE E 

k=2 1=1 p=0 q=0 



'2j'\ (2k' + 



P 



CjkC^eC^O Z{2j'-p'Z.2k'+p-qZq){2l-2.) 



(4.28) 



with j' = j — A; — 1, k' = k — i—1 and the C-symbols were defined in ( p.27| ). The action 
of Zn was defined in ( |4.25| ) (it does not act on the set of indices labeled (2£ — 2)). 



5 Curved space 

In this section we shall generalize our results to a curved space with metric g^,^. 
Riemann normal coordinates x'^ can be defined by^ 

9t^u{0) = 6/,^ , x'^gf.^ix) = x^'gi.^iO) (5.1) 



^ We use the same symbol for normal coordinates as for general coordinates. This 
should not cause confusion. In these coordinates and at the origin there is no need to 
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which is equivalent to 

9^M=^^.. , ^7Mm,/.2...M.)(0) = , J>2 (5.2) 

These properties also hold for the inverse metric g^'^ . Taking the origin of our normal 
coordinate system to coincide with the point x', we have 

a(x,0) = = \x^x''5^^ , P(a;,0) = gixY^'"^ (5.3) 



and De Witt's ansatz (|2.5| ) for the heat kernel becomes 



Kix,0;T) = (47rr)-'^/2^?(x)-i/^e-^'/^-Ea,(x,0)- (5.4) 



Using ( |5.3|) , the recursion relations ( |2.9| ) become 

{x^'d^ + j) a, = J {d^g^'d, + 2A^d^ + Z) a,_i , [ao] = I (5.5) 

Here 

Z = Z*^ + Z^ , Z*^ = X + A'^,^ + Af'A^ , Z^ = i^,/^^ - ^5,^5,^ (5.6) 
where we defined 

S = In I? = - In ^ (5.7) 

Note the position of the explicit inverse metric in (|5.5| ). The quantities Z*^ and Z^ 
are the matrix and scalar parts of Z respectively. Also note the similarity between 
- X and Z'^ . In the flat space limit Z^ vanishes and Z^'^ reduces to the quantity 
earlier defined as Z in ( [4.9[ ). We now reserve the symbol Z to mean the sum of 
and Z^. 



5.1 Scalar part 

We write Z^ with explicit inverse metric as 

Z' = lig^^'B,,),, - \ B,,g^''B,, (5.8) 

Our task is to find a covariant expression for the partial derivatives of Z^ at the 
origin of the normal coordinate system. This in turn requires the expansion for the 
(inverse) metric in normal coordinates, which is a well-known problem. We find it 
easier to work with the vielbein, the metric being defined as usual by 

g^„y = e^^ej'bah (5.9) 



distinguish co- from contravariant indices. We exploit this to position indices in such a way 
as to cause minimal clutter. Normal coordinates a;'' may alternatively be defined in terms 
of the affine connection by x^x"Y^^{x) = (which is equivalent to .,(0) = for 



j > 2). Note the similarity with the Fock-Schwinger gauge (4.4) 
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where a and h are tangent-space indices. A recursive formula for the vielbein in 
normal coordinates was given in [T^. Here we shall give its solution to all orders in 
the curvature. Using matrix notation for the vielbein, i.e. 



-/i ,111... 11 j 



(0) 



(5.10) 



and defining the following sans sen/ curvature symbol^ (compare with ( f4.12| )) 

J-1 



J + 1 



-R(mi%2'';...Mj)(0) 



J>2 



(5.11) 



we can write the recursive formula of ITOll for the vielbein as follows 
Eo 



I 



El 







k=2 



,k j(j + i: 



KfcE, , j>2 (5.12) 



where I is the dx d unit matrix and total symmetrization on the k + i = j indices is 
implied. We treat Kj as a symmetric matrix in the index pair nu. Taking the trace 
of such a matrix yields 



tr[K, 



J-1 
J + 1 



-R(/.iA»2;...Mj)(0) 



J>2 



(5.13) 



where the Ricci tensor and its covariant derivatives appear. Note that we traced over 
the world indices. This is not to be confused with the trace over field or 'bundle' 
indices which does not occur in this paper. In taking the trace over a product of 
K-matrices, say tr[KjKfcK£], total symmetrization on the j + k + i indices is implied. 
Such traces appear in our final expression for the heat kernel coefficients in curved 
space, e.g. 



tr[K(2K2)] = ii?(//(0)i?A^)^(0) = ^R^'R.x + ^R'''^''R. 
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iXfiu 



(5.14) 



where the parenthesis on the left hand side imply symmetrization and contraction 
(see ( |A.9D for the general case). In this way we maintain an index-free notation. 
Returning now to the recursion formula ( |5.12| ), the first few cases are 



■K, 



E4 = -K4 + fK(2K2) , E5 = -K5+2K(2K3)+4K(3K2) 

(5.15) 

where the parenthesis mean symmetrization only. Note that for j > 5 the Ej are not 
symmetric matrices. By iteration of ( 5.12| ) we obtain as solution to all orders 



[i/2]-i 



j-2n j2~2n+2 j„-2 



Ji = 3 



- E (-ir E E -E 

n=l ki=2 k2=2 k„=2 

3i=3 - ^ki, , i>2 



3 

.ki, .., kr, 



t\3i{3i + ^) 



(5.16) 



^°Except for the normalization factor, this agrees with the definition in [||. See (A. 2) 
for our curvature conventions. 
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where total symmetrization of the j indices is to be understood. To the best of our 
knowledge, such an explicit expression for the vielbein in normal coordinates was not 
available in the literature up to now|^. For the inverse metric it follows that 



[j72] 

g^'^j = 2 K/^ + 2 ^ 

k=2 



1 + 



ki 



jU + 1) 



(5.17) 



where £ = j — k (see ( [A.12| ) for our summation conventions). Finding the heat kernel 
coefficients through fourth order in K turns out to require knowledge of the inverse 
metric through second order only. See appendix E. Explicit expansions for the metric 
and its inverse to some finite order in normal coordinates are well known, see e.g. 
[^, ^ and our eq ( [E.4D . We next use 



B(x) 



tr \nE 



OO 1 

5:-tr[(/-i?)'=] 

k=l 



(5.18) 



Thus B and its first derivative vanish at the origin and for j > 2, after collecting 
terms of the same degree in K, we find (c/ ( [E.5|) ) 



[j72] 



\'^/ k=2 e=k 



[j/3] [(j^k)/2] 



k=2 

[i/4] [j/2]-kj-k-2i 

+ E E E 

k=2 i=k" m=k 



J ' 

.k,i,mj 



PkemntT[KkWmK] + 0(K^ 



(5.19) 



where tr is a weighted trace defined below and the lower limit k" means: when i = k, 
reduce the upper limit for the sum over m to [j/2] — k. The coefficients are given by 



Pki 


= 1 


- P2 


l^klm 


= 2 


- ^3^ 


/3k£mn 


= 2 


- P4 









k{k + l) 



2ki 



j(j + l)J j(j + l) 

k{k + l)(j + £ + l)n Akim{f + Aj + 3 + ki + im + mk) 



j{j + l)DkeDimD. 



mk 



j{j + i){j-e + i) 

k{k + 1) / n{n + + m + 1) m(m + 1) 

{k + i)Dke ^ j(j + ^)Dken 2(m + n)Dr, 

k + i + 1 , DMm = k + i + m + l (5.20) 



Here denotes the group consisting of all cyclic and anticyclic permutations on N 
objects (so IP2I =2, |P/v| = 2N for > 2) and its action is defined by 



Pjv[/(fci,...,fciv)] = E/(^- 



N 



(1)' 



(5.21) 



i=l 



where the sum runs over all cyclic and anticyclic permutations. We used here that 
the trace of a product of symmetric matrices is invariant under such permutations 



^^We have been informed by C. Schubert and U. Miiller that they have obtained results 
equivalent to ( ^.16 ). Details will be published elsewhere [ p9| . 
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in order to simplify the sums. Thus, in tr[KfcK£Km] the order of the matrices is 
irrelevant and we can assume without loss of generality that k < i < m. Similarly, 
in tr[KfcK^KmK„] we can assume k to be the smallest label and i < n (when i = k 
we arrange m < n). Ordering subtleties do not occur before fourth order, where we 
must distinguish^ between tr[K2K2K3K3] and tr[K2K3K2K3]. The weighted trace tr 
is defined by 

tr[Kfe,K,, . . . Kk^] = ^ tr[Kfc,Kfc, . . . K^^] (5.22) 

with 5" the total number of distinguishable cyclic and anticyclic permutations (includ- 
ing the identity) of the product K^^K^j . . . K^^^. We refer to appendix F for a table 
containing all symmetry factors for N < 4. Substituting ( |5.19| ) and ( |5.17| ) into 



we find for any j > but with a cutoff at third order in K (see appendix E for the 



details and ( |6.4|) for examples) 



k=2 



k=l 

ri + 2, rj-fc + 2 




-k fi)u\ 



+ E ; (7^tr[K,,M - itr[K^(,]tr[K,),: 



k + 3 




tr[K^(fc] tr[K<?Km)^] 
tr[K^(fc] K^^'^^trlK,),] 

{Kf^e^k-iVtilKj^k,.).] 



J 



K(/^tr[K,K,_fc_,^),] 



(5.23) 



This expression is complete for j < 5. The coefficient jSkim which appears here is 
obtained by replacing j by j + 2 in ( |5.20|) and {K^K^} is defined as in (|2.11| ). The 



prime on [j/2] in the second sum implies division by 2 when j is even and this limit is 
reached. Similarly, the lower limit i = k' means division by 2 when i = k. We recall 
that total symmetrization on all indices within a given trace (except for those being 
traced over) is to be understood and thus one has e.g. 

tr[K^(i]tr[K3Ki)^] = tr[K,,(i] tr[K2K2g (5.24) 



^^Thus in obtaining the logarithm of the Van Vleck-Morette determinant through ninth 
order in the normal coordinates, see (|E.5|), we may treat the Kj as commuting objects ! 
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Replacing j by 2j in ( p.23| ) and contracting fully, the second and last two terms 
vanish. Keeping also terms of fourth order in K we find (c/ (|6.5|) ) 

Z^(2i) = itr[K(2j+2)] 

+ E (^^^) (|^tr[K(,+iK,+i)] - itr[K,(,]tr[K,g 

r 2j+2 i r 2j-fc+2 i 

^ 3 I I 2 I + 2\ 1 

+ 51 XI L /) 2 

A;=2 £=fc \ ^ / 

r 2j — fc — 1 1 
2i-3 [^^] / X 

E , j7^mtr[K^(fc]tr[K™K^g 
fc=i ^=1 V^)^/ 

j-1 2j-k-2 / 2 ■ \ 

-EE /jtr[K,(,]K./^tr[K,).] 
k=i i=k' V^'^/ 

+ E E E [up 2 

„. ^ r2i-fc-2i r2j-fc-«i 

-EE E L ^ ^ h^™"t^[^^Mfc]t~r[KnK™K^g 

fc=l ^=1 m=^+l \«';<^;''V 

[¥] [2^-^] / 2^- \ 
-EE E r, ^ 27fe^7^„tr[K^(fcK£]tr[K„K,„)^] 

i-2 2i-fe-4[2^] / 27 \ 

-EE E L; a,nntr[K^(,]{K^K„r^tr[K,),] 

fc=l e=k' m=2 V^' ^' 
2i-5 [2^4^] 2j-fc-£-2 / 2,- \ 

-EE E /^)27^mtr[K^(fc]K„'^-tr[K^K,g 



fc=l ^=1 m=^+l 
'5 



+ 0(K5) (5.25) 

This expression is complete for j < 3. The a^n are defined in ( 5.17| ) (with j —* 
m + n = 2j — k — i) and the /9-coefficients are those of (|5.2CI|) with j 2j + 2. The 
7-coefficients are defined by 

- ^ - k + i + m + 1 

= k + i + 2 ' ^''^ = 2(A; + 1) ^^-^^^ 

The upper limit j' means: divide by 2 when k = j. A lower limit k' means: divide 
by 2 when i = k. The upper limit (j — k)' means: divide by 2 when i = j — k and 
m = k. Finally, the lower limit k" means: when i = k, reduce the upper limit for the 
sum over m to j — + 1. 
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5.2 Matrix part 

Abbreviate the j-th partial derivative of the gauge connection at the origin by 

A, = A^,^,...^,(0) (5.27) 

Then the recursion for the gauge connection is given by (until further notice, paren- 
thesis imply symmetrization only) 

Ao = , A, = Y, + X: (i) ^ E(,Y,) , J > 1 (5.28) 



k=2 



,k 3 + 1 



We iterate this and use ( [5.16| ) to obtain as solution to all orders 

j-2nji-2n+2 i„_i-2 / ■ \ i 1 l-. f A- ■ -^ 1 

n=l k=l ki=2 k„-i=2 V^' ••' '^n-1/ J J- j=i Jz\Jz t- IJ 

i-1 

j, = j-Y^ke , ko = k (5.29) 

where symmetrization on the j indices is to be understood. We used that, since the 
Kj are symmetric matrices, we may write them to the right of the vector Y^q in the 
reverse order. Through second order in K, ( p.29| ) reads (cf ( |E.6| )) 

A. = Y, -E('^)^Y(,K,) (5.30) 
We thus find (this is exact for j < 5; see ( p.2| ) for examples) 

k=l k=2\ J 3 + ^ 

[j/2]-lj-fc-2 / ■ \ 



fc=l £=fc 
j-2j-k 

+ EE 

fc=2 £=2 



j + l\/3£ + m + l + k{k + l] 



k,ij^i + m + l {j + 2){k + £ + l) {k + i){k + i+l)' 

x(K(fcK,)^'^Y/;,_fe_,) + 0(K3y,K2y2) (5.31) 

If we now replace j by 2j in ( p.31| ) and contract all indices, the second, fourth and 
last term vanish. Including terms of fourth order in curvatures, we find {cf (|6.3| )) 

Z^^(2i) = X(2j) + E ( / ){Y(fcY^)} 



k=l 

j-1 2j~k^2 



+ E E (:^) + JT^){v.^".Y«} (5.32) 



j_2 2i-fc-4 2i-fe-£-2 / \ mnAr« 



fc=l £=fc m=2 



,k^ i^Tfl / D Icrn.D lln.D mriD krn.ri.D P 

22 



where the D-symbols were defined in (|5.2(]| ) and 



Nkhn = 5 + 5ki + 12 j + 16 jmn + {6m + {k + m) {9m + lOn) 

+k'^{£ + m + 2n + 2) + {3k + Ai + 2m)m^ + 4Hm}s 
{/(/c, £, m, n)}s = /(fc, £, m, + /(£, /c, n, m) (5.33) 



We recognize the first fine of ( ^.32| ) as the fiat spacetime result ( |4.16|) . As we already 
mentioned there, we are not able to give a closed expression to all orders in K here. 

(2i) 



However, the terms of A^-th order in K in T}"'^ will be of the form 



N 

{Y(fcK^,...K^^Y,)} , 2j = fc + £ + ^mi , k<l (5.34) 

i=\ 

and this shows that Z^^ has a much simpler structure than Z^ . 



5.3 Expansion in derivatives 

We retrace the steps of sect 4.1, starting with (compare with ( [4.17|) ) 

aj,n = j ^ ^ ( ^i-i,^M>™ + (^^j-i),") (5.35) 

where now 

X = Z + 2A^9^ = Z^^ + Z^ + 2A''d^ (5.36) 

Except for n = or 1, we can not replace the explicit pair of indices in ( p. 351 ) by 
our shorthand notation as in ( [4.17| ). However, for even values of n and after full 
contraction we can use the following lemma: for any scalar F{x) 

F%(2n) = F(2„+2) , n>0 (5.37) 

This follows immediately from the defining properties of the inverse metric in normal 
coordinates, see (|5.2| ). Using this lemma, we obtain the following chain of substitu- 
tions (omit numerical factors and ordering labels, but keep undifferentiated Z's) 

a 3(2) + Za 

a(2n) a(2n+2) + (Z'^a)(2n) , ^ > 1 (5.38) 

Z(p^a2„-p) Z(p^a''2n-p)/i + Z(pf(Xa)2n~p) , p = 0...2n 

etc, where 

Z = Z^'' + Z^ = X + i tr[K(2)] = X(0) + ii?(0) (5.39) 

In the second step we used that ( [4.20| ), with the understanding that d"^ = d^d^, 
remains true in normal coordinates. We can remove X from the last line exactly as 
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in the first three hnes of ( 4.22|) . However, the first term on the right hand side in the 
last hne of ( ^.38| ) yields a new term as compared to fiat space {r = 2n — p) 



r + 1 
p /r + 1 



tE „ Z,(,.ith'^V._,+i). (5.40) 



^(r\ pjr-q) ^ , 
~ ~,t^oU(9 + l)(^ + 2) 

The first step merely defines h'^^ = g^*^ — 6'^'^. The vertical bars indicate that the 
index i/ is not to be included in the indicated symmetrization. In the second step we 
use that, due to (|5.2|) 

Z(;(h'^^a,|,|),^) = (5.41) 

Writing this out with respect to the position of the index /i shows the equality of the 
first two lines of ( ^.401 ). Next we use the binomial theorem and in the last step we 
shift q and use that h and its first derivative vanish at the origin. We absorb this new 
term through a redefinition of Zq in ( |4.22|) and arrive at the following conclusion. 

The diagonal values of the heat kernel coefficients in curved space are obtained from 
those in flat space by formally making the same substitutions as in l \4-24 ), where Zj 
now stands for Z^^ j + Z^j (see ( \5. 23^ , \5. 3]\ ) and the action of Z is defined by 

Z(/Z„Zfc) = Z(/Z„Zfc) + —^ZjA^n+iZk^i)^ - ^ ^^^^ ^ 2) 'Z-^^(j~h'"'n+2Zk-l)u 

(5.42) 

etc, where the (contravariant) gauge connection and inverse metric follow from l \5.11\ ) 
and (\5Jdi )- 

These redefinitions being understood, the heat kernel coefficients formally remain 
unchanged upon going to curved spacetime. 



5.4 Locally symmetric space 

The results of the previous sections nearly suffice to find 85. The few missing terms 
are of fifth order in the Yang-Mills and Riemann curvature. They are most easily 
found by considering a locally symmetric space, i.e. a space with covariantly constant 
curvatures. In that case one can find closed expressions which hold through all orders 
in these curvatures. This situation has been considered in detail by Avramidi ||24|| . 
Here, we merely want to find those terms in a given heat kernel coefficient for a general 
curved space which do not contain explicit covariant derivatives. We therefore do not 
take account of the consequences of the requirement that the Riemann curvature is 
covariantly constant (i.e. VR = would imply [V, V]-R = 0). 

24 



Thus, consider the case where only Xq, Yi and K2 are nonvanishing. Define 

Y,{x) = \f^u{Q)x^ , K^,{x) = ^/2^p..(0)x^x" (5.43) 



The vielbein is then given to all orders in K by (see (|5.16| ) 

^m-g(^(-3Ar^54^ , (5,44) 

Note that the vielbein is an even function of S and hence it depends only on K. Since 
K is a. symmetric matrix, so is E. It follows that the inverse metric simply equals 
E~'^. Differentiating i? = — tr In ii^ we find 

5,^ = i tr[ir,^L] , = I ( coth S (5.45) 

where L is an even function[^ of S. In general, the commutator of K and K^^ does 
not vanish, but the trace insures that ( ^.45| ) holds. Differentiating L we find 

L,^ = {K^^L') , {K,,K^) = - ^ K'^K^^K^-'^ (5.46) 

^ fc=0 

With this notation, the result for in a locally symmetric space reads 

Z^[K] = \{E-^Y^[ii[K^^,L+K^^L^J\-\ii[K^^L\ ii[K^,L]) + \{E--'K^^LY'' ii[K^,L\ 

(5.47) 

To obtain a similar result for we note that the covariant gauge connection is given 
by 

A[y,K\ = (^%^)'y (5.48) 

The contravariant gauge connection is therefore 

A[Y,K] = (sech 5/2)2 F ^5^49^ 
and we obtain the following result for in a locally symmetric space 

Z^[X,y,ir] = X - g-(tanh S/2f Y + y ( ^anh g/2 x 2^ ^^^^^^ 



S/2 

X + d-ijK + ^K^+ ...)Y + Y{I + ^K + + r^K^ + . . .)Y 



-'^^Except for the extra factor 3/5', this is the weU known Langevin function. 
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6 Explicit results for ai through 35 



As an application, we give here the explicit results for the diagonal values of the first 
five heat kernel coefficients, obtained from ( [4 .241 ), ( |4.28| ) and ( |5.42| ). The number of 
terms equals 1,2,4,10 and 26 respectively. The sixth coefficient in fiat space can be 
found in appendix C. 

ai = Z 

32 = + |Z(2) 

33 = Z=' + i{ZZ(2)} + iZ(itZi) + ^Z(4) 

a4 = Z^ + f{Z2Z(2)} + |zZ(2)Z + |Z(itZZi) + |{ZZ(itZi)} 

+ i{ZZ(4)} + |{Z(itZ3)} + |Z(2tZ2) + iZ(2)^ + ^Z(6) 

as = Z5 + |{Z%)} + {Z2Z(2)Z} + iZ(itz2Zi) + |{ZZ(itZZi)} + {Z2Z(itZi)} 
+|zZ(itZi)Z + f {Z2Z(4)} + f ZZ(4)Z + f {ZZ(3tZi)} + ^{Z(itZZ3)} 
+if {ZZ(itZ3)} + f Z(2tZZ2) + f {ZZ(2tZ2)} + |Z(2)ZZ(2) + |{ZZ(2)n 
+f Z^tZ^Zi) + f {Z(itZiZ2)} + |{Z(itZi)Z(2)} + ^{ZZ(6)} + ^{Z(itZ5)} 

+ ^{Z(2^Z4)} + ^Z(3^Z3) + |^{Z(2)Z(4)} + ^Z(3/Z(3) + J^^S) (6-1) 

The Z's appearing here were defined in the previous sections, but to be quite explicit 
and to avoid misunderstanding of our conventions, we give them below. We have 
Zj = Z*^j + Z^j with the matrix quantities, see (|5.31| ), given by 

Z^ = Xi+Y(/;i) 

Z^^2 = X2+Y(/,2)+2Y(iYi) + |K(2^%), 

Z^^3 = X3 + Y(/^3) + 3{Y(iY2)} + |(21K(2^^Y/;i) + 16K(3'^^Y,),) 
Z^^4 = X4 + Y(/,4)+4{Y(iY3)} + 6Y(2Y2) + 12Y(iK2Yi) 

+|(8K(2^'^Y/;2) + 9K(3'^'^Y/;i) + (5K(4 + 9K(2K2)^'^Y^),) 
Z^5 = X5 + Y(/;5)+5{Y(iY4)} + 10{Y(2Y3)} + 24Y(iK3Yi)+9{Y(iK2Y2)} 

+|(45K(2^^Y/;3) + ^K(3'^'^Y/;2) + 9 (5K(4 + 31K(2K2)'^^Y/;i) 

+4 (6K(5 + 51K(3K2 + 60K(2K3)'^'^Y^),) (6.2) 

In ( |6.2|) the parenthesis around the indices imply symmetrization only (but note that 
Z^^4 and Z^5 appear in ( |6.1| ) with at least one respectively two contracted pair(s) of 
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indices). All terms in (|6.2|) with explicit indices vanish upon full contraction. From 
(IPD we find 

Z^(2) = X(2)+2Y(iYi) 

Z^(4) = X(4)+4{Y(iY3)} + 6Y(2Y2) + 12Y(iK2Yi) 

Z^6) = X(6)+6{Y(iY5)} + 15{Y(2Y4)} + 24Y(3Y3)+40Y(iK4Yi) 

+66{Y(iK3Y2)} + 50{Y(iK2Y3)} + 72Y(2K2Y2) + 153Y(iK2K2Yi) 



1) 



Z^\s) = X(8)+8{Y(iY7)} + 28{Y(2Y6)} + 56{Y(3Y5)} + 70Y(4Y4) + 84Y(iK6Y 
+225{Y(iK5Y2)} + ^{Y(iK4Y3)} + 273{Y(iK3Y4)} + 126{Y(iK2Y5)} 

+480Y(2K4Y2) + 520{Y(2K3Y3)} + 288{Y(2K2Y4)} + if^Y(3K2Y3) 

+ 6^Y(iK3K3Yi) + 870{Y(iK2K4Yi)} + 1440{Y(iK2K3Y2)} 



5 

^f^{Y(iK3K2Y2)} + 1092{Y(iK2K2Y3)} + i:^Y(2K2K2Y2) 
+3348Y(iK2K2K2Yi) (6.3) 

The last term in "Z-^^s) was obtained from (|5.5CI| ). From (|5.23|) we find that the scalar 
quantities are given explicitly by 

Z^ = itr[Ki^^] 

Z^2 = Jtr[K2^^ + |K2KJ -itr[K,,(i]tr[Kig +3K(2'''^tr[K^g 

Z^3 = I tr[K3^^ + 4K2Ki^^] - | tr [K^(i] tr[K2g + 4K(3^^ tr[K^),] + 6K(2'^^ tr[Ki,,),] 



Z"4 = 3^tr 



7K4^^ + 40K2K2^^ + 30K3Ki^^ + 48K2K2K J - | tr[K^(2] tr[K2)^] 



-| tr[K^(i] tr[5K3)^ + 6K2Ki)^] - 6 tr[K^(i] Ks'^'^ tr [Ki),] + 10K(3'^^ tr[Ki^), 
+ (5K(4 + 36K(2K2)'^^ tr[K^),] + 2K(2''^ tr[5K2^), + GK^K^^,] 



Z^, = itr 



+ ^K2K3^^ + 15K3K2^^ + 48K2K2Ki^^ 



2 

-| tr[K^(i] tr[K4)^ + 4K2K2g - tr[K,,(2] tr[5K3)^ + GK^K,)^] 
-30tr[K^(i] K2^"^tr[K2).] - 10tr[K^(i] Ka^^^triKi),] 
+15K(2'^'^ tr[K3^), + 4K2Ki^),] + 4K(3^'^ tr [5K2^). + QK^K^).] 
+3 (5K(4 + 36K(2K2)^^tr[Ki^),] 

+6 (K(5 + 12K(2K3 + 12K(3K2)'^'^ tr[K^),] (6.4) 
where the parenthesis mean symmetrization only. Upon full contraction we find from 

(Pi 

Z^(2) = i tr[K(4) + |K(2K2)] - I tr[K,,(i] tr[Ki)^] 
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(4) 



(6) 



(8) 



14 " 



7K(6) + 40K(2K4) + 30K(3K3) + 48K(2K2K2) - f tr[K^(2] tr[K2g 



-2tr[K^,(i] tr[K3)^ + ^K^K^)^] - 6tr[K^(i] K2^'^tr[Ki), 



■tr 
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70 



140 



272, 



400, 



K(8) + f K(2K6) + ^KoKs) + ^K(4K4) + ^K(2K2K4) + ^K(2K3K3) 

+ ^K(2K2'K2)] - f tr[K^(i]tr[7K5)^ + 40K2K3)^ + 30K3K2)^ + 48K2K2Ki),, 

-f tr[K^(2] tr[K4)^ + 4K2K2)^] - | tr[5K^(3 + 6K^(iK2] tr[5K3)^ + 6K2Kig 
-3 tr[K^(i] (5K4 + 36K2K2)'^'^ tr[Ki),] - 12 tr[K^(i] K2^'' tr[5K3). + 6K2K1),] 
-60tr[K^(i] K3'^^tr[K2),] - 45tr[K^(2] K2^^tr[K2).] 

itr[K(io)] + IT^^ 



+ 11 " 



9K(2K8) + 500K(2K4K4) + 1680K(2K2'K4) + 1296K(2K2^K2) 
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12K(3K7) + 24K(4K6) + 15K(5K5) + ^K(2K2K6) + 290K(2K3K5) 
+460K(2K3K2K3)] + 540tr[K(3K3K4)] + 1788 tr[K(2K2K3K3)] 
-4 tr[K^,(i] tr[K7)^ + f K2K5);. + f K3K4);. + ^K,Ks), + ^K2K2K3);. 

+ ^K2K3K2)^ + ^K2K2K2Ki)^^ 

-14tr[K^(2] trfKg);, + f K2K4)^ + 15K3K3)^ + 48K2K2K2);.^ 
-| tr[5K^(3 + 6K^(iK2] tr[7K5)^ + 40K2K3)^ + 30K3K2)^ + 48K2K2Kig 
-f tr[K^(4 + 4K^(2K2] tr[K4)^ + 4K2K2g 

-4tr[K^(i] (7K6 + 125{K2K4} + I7OK3K3 + 720K2K2K2)'^^ trfKi^ 
-168tr[K^(i] (K5 + 12{K2K3})^'^tr[K2).] 

-f tr[K^(i] (5K4 + 36K2K2)^'^ tr[5K3), + 6K2K1),] 

-280tr[K^(i] K3^''^tr[K4). + 4K2K2).] - 42tr[K^(2] (5K4 + 36K2K2)^'^ tr[K2).] 
-24tr[K^(i] K2^^tr[7K5). + 4OK2K3). + 3OK3K2). + 48K2K2K1),] 
-112tr[K^(2] K3^'^ tr[5K3). + 6K2K1),] - 420tr[K^(2] K2'^'^ tr[K4). + 4K2K2).] 



-f tr[5K^(3 + 6K^(iK2] K2^^tr[5K3). + GK^K.^,] 



(6.5) 



where the seven terms of fifth order were found from ( ^.471 ) and also from expanding 
B through tenth order in the normal coordinates (see appendix E). 

To return to conventional notation we use 





y- \ X fr J- —T? 



(6.6) 



-(2) 
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28 



etc (these suffice for ai and 82). In 34 respectively 85 we have e.g. 



Z(itZZi) = Z;(Z(7^'^ + F^^-ii?nZ. (6.7) 
Z(2tZZ2) = iZ(2)ZZ(2) + |z.;(Z(7'^^ + 2F^^-|i?nZ.. + |z.;i?'^^'^'^Z,, 

Using such translations we have verified that our results for the first four coefficients 
agree with earlier authors, in particular and after accounting for differences 
in conventions and the occasional typographical mistake. 

Finally, to illustrate the compactness of our notation, consider the special case 
of a scalar field in a Ricci fiat space. Then only the Weyl tensor and its covariant 



derivatives can appear in the heat kernel coefficients. From appendix D of |[T^ we 
can read off the total number of general coordinate scalars in that case. Thus, for 
j = 2, 3 and 4 we expect 1, 3 and 12 terms in a^, respectively. In our notation, see 



( |5.13| ), Ricci-flatness implies that tr[Kj] vanishes. Only a few terms remain then in 



|5.23| , p.25| ). In particular, Z and Zi vanish and a2 through contain only 



32 

as 

34 



tr[K(2K2)] 

tr[K(2K4)] , tr[K(3K3)] , tr[K(2K2K2)] 

tr[K(2K6)] , tr[K(3K5)] , tr[K(4K4)] , tr[K(2K2K4)] , tr[K(2K3K3)] , (6.8) 

tr[K(2K2K2K2)] , tr[K^(iKMtr[K^Ki),] , tr[K^(iK2] triKaKi^ , ( tr[K(2K2)])' 



Here there are 1, 3 and 9 terms respectively, so that, starting with a^, our notation 
is not only index-free, but also generates less terms. This becomes even more pro- 
nounced for 35 where [1^ informs us that in the Ricci-flat case there are 67 terms. 



however in our notation there are only 17 terms. 



7 Conclusions 

We have presented the explicit diagonal values of the first five (six) heat kernel coeffi- 
cients for a general Laplace-type operator on a Riemannian (respectively fiat) space. 
To solve the pertinent recursion relations, we relied not only on well known tech- 
niques (matrix notation for field indices, Fock-Schwinger gauge and Riemann normal 
coordinates), but also used a new notation free of spacetime indices. It is this latter 
compact notation which allows us to write down the fifth and sixth coefficient in the 
first place. Insisting also on manifest hermiticity of the results, the fifth (sixth) co- 
efficient has 26 (respectively 75) terms. They were presented here for the first time. 
Beyond these coefficients, the leading terms - to the same order in derivatives or 
curvatures as needed for the fourth coefficient - for any heat kernel coefficient were 
given in the general case. To determine the heat kernel coefficients, we have found it 
useful to proceed in a few steps, starting from the simplest case of a fiat spacetime 
without a gauge field. We could show that 'turning on a gauge field' and next 'curving 



29 



spacetime' is taken care of by specific covariant substitutions in the fiat space coefii- 
cients. With this 'dressing up' understood, the number of terms and their numerical 
prefactors do not change in the process. 

In our notation, a typical term in a heat kernel coefficient consists of a maximally 
symmetrized (product of) covariant derivatives of the basic curvatures, made into a 
scalar by contracting all indices. Consider e.g. the following term 

tr[K(3K3K4)] = ^R''i/.;xRW;,RW;pp) (7.1) 

which appears in the fifth coefficient. Since we do not integrate over spacetime, 
partial integrations are not permitted. Furthermore, as long as we do not write out 
the symmetrization, the Bianchi identity can not be used here. Thus our notation 
gives a certain degree of uniqueness to the appearance of the heat kernel coefficients, 
missing in a more conventional notation with explicit spacetime indices. 

Of course, depending on the application one has in mind, our notation may or 
may not be useful. We plan to calculate the chiral anomaly based on our results, i.e. 
essentially evaluate the spinor-trace tr[72j+iaj] in d = 2j dimensions. In this case 



elegant and complete results are known ||2^, using the language of differential forms. 



It should be interesting to see how this can be related to our use of symmetric tensors. 



Also the gravitational anomalies p6| in = 4j ' + 2 are computable in our framework. 
Possibly, the absence of certain terms from the heat kernel coefficients found here, 
see ( p.l9|) , plays a role in this connection. 



Another task, now in progress is to work out the functional trace of the 
diagonal heat kernel coefficients given here. In that case we can maintain our index- 
free notation (e.g. the integral of ( [7.1|) is easily seen to vanish without writing out 
the indices). Although it would probably hold few surprises, c/ 1|18[, the result for the 
integrated and traced fifth coefficient could be used to study for the first time the 
one-loop short distance divergences of ten-dimensional supergravity. 



As this paper neared completion, we were informed by the authors of ||TT| that 
they had extended their results so as to include gauge fields. In they present the 
functional trace of the first six heat kernel coefficients for this case in fiat spacetime, 
reduced to a so-called minimal basis. It should be possible to compare their result to 
the sixth heat kernel coefficient as presented here, after taking its trace and integrating 
it. 

Finally, to avoid misunderstanding, we should mention that a recent preprint 
with the title "The heat kernel coefficient on a manifold with boundary" is 
not concerned with the fifth heat kernel coefficient as presented here. Rather, on a 
manifold with boundary the expansion ( p.5|) is in powers of a/t rather than r and 
therefore our aj corresponds to a2j of [^. Thus, in the terminology of we have 
determined here 'the volume part of aio'. 

I would like to thank J.-P. Bornsen for his assistance in using Mathematica to 
verify some of the results presented here. I am also grateful to I. Avramidi for useful 
correspondence. 
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A Notation 



We use Greek letters for spacetime indices. Indices j, k, £, . . . simply enumerate 
various objects. We take differential operators to act on everything to their right. 
If this is not intended, we use comma (semicolon) notation for partial (covariant) 
derivatives 

F, ^ [d,,F] , F.,, ^ [V„F] (A.l) 

A semicolon denotes simultaneous gauge and gravitational covariant differentiation. 
We use a = 6 to define a in terms of b. Our curvature conventions are fixed by (Va is 
a gauge singlet) 

[Viy,V^]VA = 'iVx-l^u] = R'^X^ivVk , R\v = R^Xkv , R = -R'^A (A. 2) 

In Fock-Schwinger gauge and Riemannian coordinates, see ( [4.4| ) and ( [5.1[ ), we haveQ 

F^,...^,(0)=F;(^,...^,)(0) , (A.3) 

for any, possibly matrix valued, general coordinate scalar F{x). This property allows 
one to immediately covariantize the partial derivatives. Due to the total symmetry 
here it is convenient and sufficient to keep track of only the number of indices, leading 
us to define the sans sen/ symbols 

F, = F;(;.....M,)(0) (A.4) 

We thus use an index-free notation for totally symmetrized tensors. At the same time, 
such a symbol implies evaluation at the origin. Following Avramidi, we use parenthe- 
sis around an even enumerative label to indicate not only total symmetrization, but 
also full contraction as in 

F(2i) = F,i,,,,...,,,,){0) (A.5) 

Since the metric at the origin is flat there is no need to write the indices in their 
covariant respectively contravariant positions. Par ahus de language an odd label 
enclosed in parenthesis denotes symmetrization and contraction of all but one of the 
involved indices. This will occur only if there is exactly one other factor with an odd 
label, as e.g. in 

F(3)G(3) = F;(.m(0)G;(«,,)(0) (A.6) 

We generalize Avramidi's notation by allowing such a simultaneous symmetrization 
and contraction to extend over several factors, e.g. 

F(.G2,_.) = F;(,,...,,(0)G;,,,,...,,^,)(0)5^^^^ . . . 5^-^'^>^-^ (A.7) 

If desired, such an expression can be written out in such a way that the total sym- 
metrizations extend only over the individual factors, as e.g. in 

F(2G2) = |f^(^;.)(0)G;(..)(0) + |F;(^,)(0)G;(,.)(0) (A.S) 

^"^Parenthesis around j explicit indices denotes total symmetrization (with division by 
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In general, expanding (|A.7|) in this way yields at first (2j — 1)!! terms. However, due 
to the total symmetry of the covariant derivatives acting on each factor and assuming 
that k < j, there remain only [/c/2] + 1 terms, [k/2] being the maximal number of self- 
contractions possible for F^. Finding the coefficient of each term is a combinatorical 
problem with the following solution 



F(2jG2fc) 



(2j + 2k 



1 V P^^' P2fc (nf\\ c{2j-2l) r(2k-2e) 

,^ _ Z^-^i-r ^tii...ti2t %i...A'2^ 



i=0 

1 ^ 



PI = (si) (2^-1)" ' ^^[^'/S] , (-1)!! = 1 (A.9) 

Here Pi is the number of ways in which one can choose k pairs out of j objects. The 
covariant derivatives on F respectively G on the right hand side are understood to 
have been totally symmetrized. In general there will be more than two factors, but 
in practice these are easily taken care of, as e.g. in 

12 8 
F(2G2H2) = Y^F(2)G(2)H(2) + Y^(F(2)G;(^,,)H.(^j,) + 2 more) + Y^F.(^i,)G;(;,p)H.(p^) (A.IO) 

Multinomial coefficients are defined as usual by 

L \ )= TtA. f ' ^^+1 ^J-T.kn (A.ll) 

\ki,...,kNj ki\...kN+il 

If the index /cat+i appears in the summand of a sum involving this multinomial co- 
efficient, then it is understood to have the indicated value. This convention applies 
also to alphabetic ordering as e.g. in 



T.( JF^Ge ^Ei JFkG,., (A.12) 

fc=0 \ / A;=0 \ / 



B Moving boxes 



In the main text we gave an expression, eq (|3.31| ), for the diagonal values of the heat 
kernel coefficients aj through third order in the matrix potential X which holds only 



for j < 5. We indicated in ( p.30|) how to obtain ae as well without generating terms 



with overlapping derivatives. In the general case we rely on the following lemma to 
move the boxes from G to F 

n 

F(2nG2p){2q-2p) = ^M^P^ (l)P{2n-2j){2jG2q) , n<p<q (B.l) 

j=0 



/3,(n,p,g) = , /(j,p,g) = (2j + 2g-l)!!P,^!,. 

* n-j) f{n,q,p) 
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with a similar expression for the case -F(2„+iG'2p+i)(2g-2p-2)- Note that this can only 
be done when all indices are contracted. To prove the lemma, expand both sides 
using ( |A.9| ) and compare. We also note the following exceptional case 

F{k+2nGk){2n) = F{2n){kG k+2n) (B-2) 



C 36 in flat space 



Expressions ( p.20| ) , (|3.31|) , ( |4.24| ) and ( [4.28| ) suffice to write down the sixth heat kernel 
coefficient in a flat space, but with a gauge connection. We obtain 

aa = + f {Z^Z(2)} + f {Z^Z(2)Z} + ^^Z^Z,,,Z^ + f Z(,tZ^Z,) + |{ZZ(/Z^Z,)} 
+f {Z^Z(itzZi)} + f ZZ(itZZi)Z + liZ^^i)} + f{Z'Z^ni)Z} + ^{Z^Z(4)} 
+^{Z^Z(4)Z} + f {ZZ(itZ3)Z} + if {Z^Z^tZs)} + f {Z^Z^tZi)} + f {ZZ(itZZ3)} 
+f {ZZotZZy} + ^{Z(itZ^Z3)} + ^{Z%n2)} + f {Z^Z(2)n + ^2(2^2^22) 
+f Z(2)Z^Z(2) + ^{ZZ(2tZZ2)} + |{ZZ(2)ZZ(2)} + f5zZ(2tZ2)Z + |ZZ(2)^Z 
+^{Z(itZiZZ2)} + |{Z(itZi)ZZ(2)} + ^{ZZ(itZiZ2)} + f {ZZ(itZi)Z(2)} 
+g{ZZ(2^ZiZi)} + |{ZZ(2)Z(itZi)} + f {Z(/ZZiZ2)} + f {Z(itZZi)Z(2)} 
+^{ZZ(itZ2Zi)} + |{Z(itZZ2Zi)} + ||Z(itZiZiZi) + |Z(itZi) Z(itZi) 
+ A{Z2Z(6)} + ^ZZ(6)Z + |-{Z(itZZ5)} + ^{ZZ(itZ5)} + ^{ZZ(5tZi)} 

+ ||{Z(2^ZZ4)} + f {ZZ(2^Z4)} + f {ZZ(4^Z2)} + f {ZZ(3^Z3)} + |lZ(3^ZZ3) 

+ |{Z(;ZiZ4)} + f Z(itZ4Zi) + ^{Z(itZ2Z3)} + |^{Z(itZ3Z2)} + f {Z(2tZiZ3)} 

+ ^Z(2'^Z2Z2) + ^{ZZ(4)Z(2)} + ^{Z(2)ZZ(4)} + ^{ZZ(2)Z(4)} + |{ZZ(3/Z(3)} 

+ f Z(3/ZZ(3) + |{Z(,tZi)Z(4)} + ^{Z(itZ2)Z(3)} + f {Z(2tZi)Z(3)} + ^{Z(2)Z(/Z3)} 

+ f {Z(itZ3)Z(2)} + f {Z(2Z2)Z(2)} + iZ(2)3 + ^{ZZ(8)} + §^{Z^,%)} + i{Z(2^Z6)} 

+ ^{Z(3^Z5)} + |Z(4"^Z4) + j|g{Z(2)Z(6)} + ^{Z(3/Z(5)} + ^Z(2)(2'^Z2)(2) + 4^Z(io) 

(C.l) 



The Z's appearing here were defined in ( |4.15| ), ( [4.16| ) and ( |4.25|) . Note in particular 
that Z = X. There is a total of 75 terms. In curved space, ag will look exactly the 
same, but in that case we do not know the values for some of the Z's. 
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D Consequences of Fock-Schwinger gauge 



To prove ( [4.20|) , we start from 

(Am^),M2...^.)(0) = (D.l) 



which holds for any function F in Fock-Schwinger gauge, see (4^)- Now take j = 
2n + 2, contract all indices and write out with respect to the index of the gauge 
connection. It is essential for the proof that, due to the full contraction, all indices 
here are dummies. 

A further consequence of ( [4.5| ) is the following. If we define Z to act to its left as 
follows 

Z(;Z„t = Z(;Z„t - JL-Z.^.jYn+i" (D.2) 

where we left out inessential factors to the right, then we have the following alternative 
notations 

Z(/Zj"Zfc) = Z(/Z„Zfc) 
Z(/Z;Z„Z,) = Z(;ZpZ„Z,) (D.3) 

In the last case the 'check' notation has the advantage that it generates less terms 
than the 'hat' notation, namely 

- " 2j 2k 

'^i}^v^-J-k) = Z(/Zp''Z„Zfc) — Z^(j„/Yp+i^Z„Zfc) + ^ ^ ^ Z(^^Z^\ n+i'^k-iy 

■ Z^(j-i^Yp+i^Y„+i'^Zfc_i)i, (D.4) 



(p+ + 1) 



E How to get your Z's 



Here we shall give some details on how to obtain ( ^.23| ). The steps for ( |5.31|) are very 



similar so we omit them. Our starting point is (|5.8| ). We take j partial derivatives of 
this expression and evaluate at the origin {m = j — k — i) to obtain 

J /j _|_ A j-ij-k / j \ 

Z^j = 2 H ( ^ jg'"'{kQj-k^t)u - L A^t^{kg'"'mQi)u (E.l) 

where we used that B and its first derivative vanish there. Separate off the terms 
with undifferentiated inverse metric 



Z^i - l^t^t^j + IJ2 ^ ^g^''ik^j-kt,)u 



i-i /j\ i-3j-fc-2 / j \ 

"4 II (^^j^MfcB^V - 4 II J2 ^j^t^{kg'"'m^e)u (E.2) 
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Inserting ( 5.17| ) and ( 5.19 ) and collecting terms of the same order in the curvatures 
K and Y, we obtain ( |5.23|) . Up to here all contractions were explicitly indicated and 
the parenthesis denoted symmetrization only. If we now replace j by 2j and fully 



contract, the second term in ( [E.2|) vanishes due to (|5.2|) and we find 

2j-3 2j-k-2 



(2i) 



k=l \ ^ / " k=l e=i 



IE 



2j 



This yields ( ^.25 ). Thus we see that, since B is at least of first order in curvature, 
finding the heat kernel coefficients through order n in the curvature requires the 
inverse metric only through order n — 2. 

Writing G for the inverse metric, we find through sixth order in the normal coor- 
dinates 



G{x) 



1+ 2K 



X 

'21 



2K 



X 

'3[ 



2(K4 + fK2^)^ + 2(K5 + 12{K2K3})|^ 



.6 



+2(K6 + if {K2K4} + ^Ks' + ^K/)|- + Oix 



7\ 



(E.4) 



For the logarithm of the Van Vleck-Morette determinant we find through tenth order 
in the normal coordinates 



B{x) 



tr 



9 d ^ 



+ (K6 + ^K2K4 + f K3' + #K2^)?7 + (K7+i^K2K5 + 15K3K4 + 48K2'K3^'' 



7 



28, 

3 

56, 

5 



5 ^9! 



70, 

3 

168, 

5 
9 



Ik 

9 



7! 



+ (K8 + 4?K2K6 + ^K3K5 + if K4^ + ^K2^K4 + i?^K2K3^ + ipK2^)^ 



+(K9 + ^K2K7 + ^KaKe + 56K4K5 + '^-^K2^K^ + 584K2K3K4 + I44K3' 



+(K,o + ^K2K8 + fiK3K, + »K,Ke + ^^k,^ + l^K2^Ke + i^K2K3K, 
, 8000, 



11 ' 11"'^"' ' 11 -'^-'> ' 11 ' 55 "° ' 11 

K2K4' + IO8OK32K4 + ^^K2'K4 + 3576K2'K32 + 1^K2K3K2K3 



11 

20736(^^5 



11 



11 



'To! 



+ 0{x'') 



(E.5) 



For the gauge connection (with covariant index) it suffices to know that 

„2 o 3 10 8 

A{x) = Yix^ + Y2- + (Y3 - 2 K2Yi)^ + (Y4 - -K2Y2 - 5 K3Yi)^ 

+ (Y5 - YK2Y3 - 5K3Y2 - |(5K4 - 9K2')Yi)^ 



3 

+ (Y6 - y K2Y4 - y K3Y3 - ^(5K4 - 9K2')Y2 



12 



-(K5-2K3K2-4K2K3)Yij- + 0{x 



X 



(E.6) 



To convert the above expansions to conventional notation, use ( [4.12 ) and ( ^.11| ). 
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F Symmetry factors 



Below we give a table containing all symmetry factors S/Pn, defined in ( |5.22|) , for 
< 4. The list ki . . . in the second column is an abbreviation for tr[Kfcj . . . K^^^]. 
It is to be understood that k, i, m and n are different integers, each having at least 
the value 2. 



N 


ki . . . k^ 


S 


S/Pn 


2 


k k 


1 


1/2 




k e 


2 


1 




k k k 


1 


1/6 


3 


kk£ 


3 


1/2 




k i m 


6 


1 




k k k k 


1 


1/8 




kiki 


2 


1/4 




kkki 


4 


1/2 


4 


kkii 


4 


1/2 




k i k m 


4 


1/2 




k k i m 


8 


1 




k i m n 


8 


1 
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